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Abstract 

We present high-resolution spectroscopy of gaseous CO absorption in the fundamental ro-vibrational 
band toward the heavily obscured active galactic nucleus (AGN) IRAS 08572+3915. We have detected 
absorption lines up to highly excited rotational levels (J ^ 17). The velocity profiles reveal three distinct 
components, the strongest and broadest (Av > 200 km s" 1 ) of which is due to blucshiftcd (—160 km s" 1 ) 
gas at a temperature of ~ 270 K absorbing at velocities as high as —400 km s _1 . A much weaker but even 
warmer (~ 700 K) component, which is highly rcdshiftcd (+100 km s _1 ), is also detected, in addition to 
a cold (~ 20 K) component centered at the systemic velocity of the galaxy. On the assumption of local 
thermodynamic equilibrium, the column density of CO in the 270 K component is Nqo ~ 4.5 x 10 18 cm" 2 , 
which in fully molecular gas corresponds to a H2 column density of A% 2 ~ 2.5 x 10 22 cm -2 . The thermal 
excitation of CO up to the observed high rotational levels requires a density greater than n c (H 2 ) > 2 x 10 7 
cm" 3 , implying that the thickness of the warm absorbing layer is extremely small (Ad < 4 x 10~ 2 pc) even 
if it is highly clumped. The large column densities and high radial velocities associated with these warm 
components, as well as their temperatures, indicate that they originate in molecular clouds near the central 
engine of the AGN. 
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1. Introduction 

Active galactic nuclei (AGNs) are broadly classified into 
two types: type 1 AGNs, which display broad permitted 
emission lines in the optical spectrum; and type 2 AGNs, 
which display both permitted and forbidden lines with 
narrow line widths. The detection of broad emission lines 
in polarized light from type 2 AGNs (Antonucci & Miller 
1985) present a clue to tie these two types together. The 
AGN activity is powered by a supermassive black hole 
(10 6 -10 9 Mq) and its accretion disk, which extends to 
<1 pc. This central engine is surrounded by an optically 
and geometrically thick dusty torus, extending to ~ 1- 
100 pc. In the "unified model" of AGNs (e.g., Antonucci 
1993), it is postulated that much of the variety in types 
of AGNs is not due to intrinsic differences between them, 
but is largely the result of their strongly non-spherical 
geometries and the different orientations of their central 
regions with respect to our lines of sight. Specifically, in 
a type 1 AGN the torus is observed face-on and the inte- 
rior may be observed directly, while in a type 2 AGN the 

* Based on data collected at Subaru Telescope, which is operated 
by the National Astronomical Observatory of Japan. 



torus is observed edge-on, which prevents a direct view 
of the nucleus. The obscuring torus around the central 
engine is the key element for the AGN unified model, and 
the strongest verification of it would come from the de- 
tection of the molecular torus itself. Recent observations 
at many wavelengths, for example the detection of silicate 
emission from type 1 AGNs (e.g., Hao et al., 2007), surely 
show the presence of obscuring tori surrounding the nu- 
clei. However, the physical conditions and structure of 
the torus have never been measured directly, and we have 
only a crude idea of them. The direct measurements of 
the physical conditions (temperature and column density) 
of the molecular torus and the determination of the gas 
kinematics within the torus are long awaited goals of ex- 
tragalactic astronomy. 

Carbon monoxide (CO) is the most abundant inter- 
stellar molecule after hydrogen molecule (H2). For the 
study of molecular clouds, pure rotational emission lines 
of CO have been extensively observed in the millimeter 
and sub-millimeter wavelength regions. However, to date 
such observations have not been suitable for the detailed 
characterization of distant molecular clouds around the 
AGNs, because they have provided information on neither 
the detailed spatial structure (due to the relatively large 
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beam sizes of current millimeter telescopes) nor the phys- 
ical conditions (due to the limited number of observable 
lines) of the clouds. In addition, millimeter observations 
can be contaminated by large-scale CO emission in the 
host galaxies on scales of 10-100 kpc, which can make it 
difficult to clearly distinguish the molecular gas in close 
proximity to the AGN, i.e. the putative molecular torus. 

The technique we employ in the current work is high- 
resolution spectroscopy of absorption lines of the CO fun- 
damental (v=l<— 0) ro-vibration band, whose center lies 
near 4.7 /im. Continuum emission associated with the 
bright, compact central region of an AGN is used as a 
background source and the foreground molecular clouds 
are to be observed in absorption. This technique achieves 
very high spatial resolution, since the resolution is deter- 
mined by the size of the background continuum source. In 
addition, a large number of CO lines covering a wide range 
of rotational levels can be observed simultaneously with 
the same instrument under the same conditions. This en- 
ables accurate determination of the physical conditions, 
such as temperatures and column densities of the absorb- 
ing molecular clouds. 

In spite of these advantages, most previous ground- 
based observations of the 4.7 /im CO absorption spec- 
trum have been limited to Galactic sources (e.g., Scoville 
et al. 1983; Mitchell & Maillard 1993; Nakagawa et al. 
1997; Goto et al. 2003), with only a few examples toward 
extragalactic objects (Spoon et al. 2003; Geballe et al. 
2006). This is mainly due to difficulties in achieving high 
sensitivities at the wavelength around 4.7 /im. 8-m class 
telescopes with high sensitivity, now are making it possi- 
ble to observe the CO ro-vibrational absorption lines in 
molecular clouds around some AGNs. 

One class of AGNs whose central nuclei are obscured by 
gas and dust in our line of sight, the so-called Seyfert 2 
galaxies, would appear to be the most suitable for such ob- 
servations. However, spectra of objects in this class have 
to date shown no significant absorption by the CO funda- 
mental ro-vibration band (Lutz et al. 2004), in some cases 
down to very low limits (Mason et al. 2006; Geballe et 
al. 2009). On the other hand, the Spitzer Space Telescope 
has detected a strong absorption by gaseous CO toward 
the dusty Ultra-Luminous InfraRed Galaxy (ULIRG), 
IRAS F00183-7111 (Spoon et al. 2004). Observations 
at the VLT by Risaliti et al. (2006) and Sani et al. (2008) 
also have detected CO fundamental band absorption to- 
ward some ULIRGs whose central sources are heavily ob- 
scured. These results suggest that CO absorption does 
not appear in typical Seyfert 2 galaxies but can be found 
in at least some heavily obscured AGNs. The reasons for 
this are unclear, although in the case of the prototypical 
Seyfert 2 galaxy, NGC 1068, Geballe et al. (2009) have 
concluded that most of the gas in front of its AGN is 
diffuse in nature and therefore contains little CO. Both 
classes of objects, Seyfert 2 galaxies and ULIRGs, ap- 
pear to have considerable obscuring material near their 
AGNs. However, hard X-ray observations (Franccschini 
et al. 2003; Ptak et al. 2003) suggested that AGNs in 
ULIRGs are not necessarily the dominant source of nu- 



clear luminosity and that nuclear star formation is a sig- 
nificant contributor. Ongoing star formation in the nuclei 
of ULIRGs implies the presence of large amounts of molec- 
ular gas there. 

ULIRGs (Sanders & Mirabel 1996) are galaxies whose 
extremely large, quasar-like luminosities (Lir > 10 12 I/ Q ) 
emerge at least 90 percent in the form of infrared dust 
emission. They are associated with interacting or merg- 
ing gas-rich galaxies (Armus et al. 1987; Sanders et al. 
1988; Leech et al. 1994; Clements et al. 1996; Murphy 
et al. 1996), in which the collision produces either an ex- 
treme burst of star formation (usually near the nucleus of 
one of the merging galaxies), a greatly increased rate of 
infall into the vicinity of an AGN, or both. In many cases 
it is difficult to distinguish between the two phenomena, 
because the large amount of dust prevents direct observa- 
tion of the energy source. However, ULIRGs with bright 
pointlike nuclei at mid-infrared wavelengths, whose spec- 
tra show strong dust absorption, are generally thought 
to possess luminous buried AGNs (Imanishi, Dudley & 
Maloney 2006). In the case of a ULIRG-producing merger 
involving a galaxy that contains an AGN, the survival of 
the molecular torus in the AGN is unclear. The torus 
could be disrupted by gas falling into the nucleus or by 
ejection of gas caused by the rapid increase in luminos- 
ity of the AGN precipitated by such infall, or simply by 
intense heating by the AGN alone. 

Infrared spectroscopy of ULIRGs is a possible way to 
investigate the energetic events in the central region of ob- 
scured AGNs. For such purpose, the observations of the 
CO fundamental band absorption toward ULIRGs with 
space telescopes, Spitzer and AKARI, are effective in pro- 
viding some information (Spoon et al. 2005, Imanishi et 
al. 2008, Imanishi et al. 2010). However, the spectral res- 
olutions of those facilities are relatively low (R ~ 600 for 
the Spitzer/IRS, R ~ 100 for the AKARI/IRC) and are 
insufficient to resolve molecular ro-vibrational bands into 
individual transitions, let alone resolve velocity structure. 
Thus, in order to reveal the physical conditions of the cir- 
cumnuclcar region and the origin of the CO absorption, 
we have obtained high-resolution spectra toward heavily 
obscured AGNs, using ground-based telescopes. 

In this paper we report on the high resolution 
4.90-5.13 /im spectrum toward the obscured AGN 
IRAS 08572+3915 NW, using the 8.2-m Subaru Telescope. 
The detection of CO absorption lines in the fundamen- 
tal band toward this galaxy was reported by Geballe et 
al. (2006), using the United Kingdom 3. 8-m Infrared 
Telescope (UKIRT) . Although the signal-to- noise ratio of 
their spectrum is modest, they detected broad and com- 
plex absorption line profiles. Their spectrum shows that 
the line profile are dominated by a very broad absorption 
feature, which extends from to —400 km s _1 and peaking 
near —160 km s _1 relative to the systemic velocity. They 
also detected a second velocity component close to the 
systemic velocity, and marginally detected a third compo- 
nent that is rcdshiftcd. Our new and improved spectrum 
covers a greater portion of the band and has a consider- 
ably higher signal-to-noise ratio. Using it we can clearly 
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separate the CO absorption into the three components de- 
duced from the earlier spectrum. The new spectrum can 
be used to better constrain the physical conditions in the 
absorbing molecular clouds. 

The characteristics of IRAS 08572+3915 are reviewed 
in the next section. The observations and the data reduc- 
tion procedures are described in Section 3. In Section 4, 
we present the observed spectrum and the line profiles. 
Temperatures and column densities of the absorbing gas 
are estimated in Section 5.1 assuming local thermody- 
namic equilibrium (LTE) . The origin of each velocity com- 
ponent of the absorbing gas is discussed in Section 5.2, 
and the line of sight dimensions of the warm molecular 
clouds is estimated in Section 5.3. The locations of the 
continuum emitting and absorbing regions are discussed 
in Section 5.4. In Section 6, we discuss the relation be- 
tween the warm molecular clouds we discovered and the 
putative AGN molecular torus. A summary is given in 
Section 7. 

2. IRAS 08572+3915 

IRAS 08572+3915 is a warm ULIRG (im( 8 -iooo /1 m) = 
1.21 x 10 12 L Q , /25//6O = 0.23), discovered by IRAS 
(Sanders et al. 1988). The redshift of this galaxy, deter- 
mined from CO line emission at millimeter wavelengths is 
0.0583 (Evans et al. 2002). The galaxy is optically clas- 
sified as a Low-Ionization Nuclear Emission-line Region 
(LINER) galaxy by Veilleux et al. (1999). The Hubble 
Space Telescope (HST) NICMOS image shows double nu- 
clei separated by 5.5" (~ 6.6 kpc), with the northwest 
nucleus overwhelmingly the more luminous (Scoville et 
al. 2000). Mid-infrared high- resolution images indicate 
that almost all the mid-infrared luminosity comes from a 
spatially highly compact source (< 0"3) in the northwest 
nucleus (Soifer et al. 2000). The mid-infrared 8-13 /xm 
spectrum shows an extremely deep silicate absorption fea- 
ture (T9.7 ~ 5.2; Dudley & Wynn- Williams 1997, r 9 . 7 ~ 
4.2; Spoon et al. 2006, Spoon et al. 2007), and the near- 
infrared 3-4 lira spectrum is dominated by the 3.4 /xm 
hydrocarbon absorption feature (T3.4 ~ 0.9) with no ev- 
idence for emission by polycyclic aromatic hydrocarbons 
(PAHs) (Imanishi & Dudley 2000; Mason et al. 2004; 
Imanishi, Dudley & Maloney 2006; Geballe et al. 2009). 
On the basis of their spectra, Imanishi, Dudley & Maloney 
(2006) argued that the energy source of IRAS 08572+3915 
is very compact and is surrounded by optically thick dust, 
as is expected for a buried AGN. On the other hand, 
IRAS 08572+3915 does not show any indication of AGN 
activity in X-rays (Risaliti et al. 2000) or at radio wave- 
lengths (Smith, Lonsdale & Lonsdale 1998). This implies 
that the AGN is deeply buried in gas and dust (iVn ^> 10 24 
cm" 2 ), i.e., it is a Compton thick AGN. 

3. Observations and data reduction 

M-band spectra of the northwest nucleus (NW) of 
IRAS 08572+3915 were obtained using the Infrared 
Camera and Spectrograph (IRCS; Tokunaga et al. 1998; 
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Fig. 1. Pair-subtracted M-band image of IRAS 08572+3915 
with position of slit indicated. The telescope was moved to 
alternately put the source at the A and B positions in order 
to facilitate subtraction of the sky emission. 

Kobayashi et al. 2000) on the 8.2 m Subaru Telescope 
atop Mauna Kea, Hawaii. The observations were per- 
formed on January 15 UT, 2004 for the 5.00-5.13 /im 
spectrum and on February 3 UT, 2004 for the 4.90-5.04 
/zm spectrum. Seeing conditions were 0"6— 0"8 FWHM 
in the R and / bands, in which we monitored the see- 
ing, for both nights. An echelle and a cross-dispersing 
grating in the IRCS were employed to cover the following 
lines of the 12 C 16 fundamental band in one exposure: 
P(7)-P(19) on January 15, R(3)-R(0), and P(l)-P(ll) 
on February 3. A slit width of 0"6 corresponding to a 
resolving power of R ~ 5,000 (Aw ~ 60 km s _1 ) was em- 
ployed. The slit length was 10"5. The slit position angle 
(P. A.) was set to 55° east of north in order to observe only 
the northwest nucleus and avoid contamination from the 
southeast one. Figure 1 shows the slit position overlaid 
on the M-band image obtained with the slit viewer cam- 
era. The telescope was nodded 5."0 along the slit with the 
nodding frequency of ~ 0.015 Hz for subtraction of the 
sky emission and the dark current. The on-source inte- 
gration times were 64 and 44 minutes on January 15 and 
February 3, respectively. 

In order to cancel the telluric absorption, we obtained 
spectra of the early type standard stars HR 4534 (A3 V; 
V=2.14) and HR 3982 (B7 V; V=l.35) immediately after 
each observation with an airmass difference less than 0.15. 
We also obtained a spectroscopic flat field at the end of 
each night using the halogen lamp in the calibration unit 
in front of the entrance window of IRCS. 

Standard data reduction procedures were employed. 
The observed spectral frames in the nodding sequence 
were coadded and differenced in order to remove the at- 
mospheric emission. Flat-fielding was performed by di- 
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Fig. 2. Spectrum of IRAS 08572+3915 NW showing fundamental band ro-vibrational lines of CO at 4.900-5.128 fim. Data points 
shown in gray are excluded from the analysis, because of the low atmospheric transmission (< 0.5%) or low signal to noise ratio (> 
10%). The adopted continuum is drawn as a thick line. The wavelengths of the R- and P-branch lines of 12 C 16 corresponding to 
the redshift of the host galaxy are indicated on the top by solid vertical lines. The error is estimated from the difference between the 
spectra obtained at the A and B positions. The atmospheric transmission curve was calculated with the ATRAN code (Lord 1992). 



viding by the dark-subtracted lamp spectral image. Bad- 
pixels were identified on the basis of the statistics of the 
dark and flat field images, and removed by interpolation 
over surrounding pixels. One-dimensional spectra were 
extracted using the aperture extraction package in IRAF 1 . 
The emission from IRAS 08572+3915 NW was unresolved 
along the slit. The widths of the extraction apertures were 
adjusted to maximize the signal-to- noise ratios of the spec- 
tra and were 0"70 on January 15 and 0"46 on February 

3. Wavelength calibration was performed by comparing 
the positions of atmospheric absorption lines in observed 
spectra with the wavelengths of the lines in a model atmo- 
spheric transmission curve calculated using the ATRAN 
code (Lord 1992). The accuracy of the wavelength cali- 
bration is better than 6 km s . 

4. Results 

The observed M-band spectrum of IRAS 08572+3915 
NW is shown in Figure 2 together with the ATRAN model 
spectrum of the atmospheric transmission. Absorption 
lines of CO in the v=0 state arising from rotational levels 
as high as J=17 were detected. These lines are extremely 
broad; the typical full widths at zero intensity (FWZIs) 
are more than 500 km s" 1 . The maximum depth of the 
absorption lines is about 60% of the continuum level. This 
implies that the area covering factor (C/; see Figure 3), 
the fraction of absorbed area in the line of sight, is at least 
60%. 
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IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatories (NOAO) in 
Tucson, Arizona. NOAO is operated by the Association of 
Universities for Research in Astronomy (AURA), Inc., under 
cooperative agreement with the National Science Foundation. 
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Fig. 3. Schematic view of a covering factor (Cf), the frac- 
tion of the continuum source obscured by absorbing clouds 
along the line of sight. (Left) covering factor is unity. (Right) 
covering factor is 0.6. 

Calculated wavelengths of the R- and P-branch lines 
of the w=l+-0 band of 12 C 16 0, shifted to the redshift of 
the host galaxy, are denoted in Figure 2 by solid vertical 
lines. The spacings of the observed lines agree well with 
the calculated ones. This wavelength region also contains 
absorption lines of isotopomers, 13 C 16 0, 12 C ls O, etc., and 
of other vibrational transitions such as v—2<— 1, w=3<— 2, 
and so on. There is no indication that any of those bands 
are present. The strengths of the isotopic lines would be 
much less than those of the w=l+-0 lines of 12 C 16 for the 
standard isotope ratios (Sakamoto et al. 1997) unless the 
12q16q i mes are extremely heavily saturated. The spac- 
ings of lines of other isotopomers and other vibrational 
states of 12 C 16 do not match those of the observed lines. 
Therefore, we conclude that all the observed absorption 
lines are w=l+-0 lines of 12 C 16 0. In the following, we use 
CO for 12 C 16 0. 

Figure 4, which shows the observed velocity profiles for 
various transitions, reveals remarkably broad and high- 
speed components. The velocity scale takes into account 
the orbital motion of the Earth relative to the Sun and 
the motion of the Sun relative to local standard of rest 
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Fig. 4. Observed line profiles for various transitions in the 
spectrum shown in Figure 2. The horizontal axis is the veloc- 
ity relative to the systemic velocity of the host galaxy. 

(Meeks et al. 1976). Although the line profiles are com- 
plex, they can be generally described as consisting of three 
velocity components: one centered at km s _1 (the sys- 
temic velocity); and the others peaking at —160 km s" 1 
and +100 km s _1 . The relative intensities, velocity ex- 
tents, and degrees of overlap of these components with 
one another vary with rotational level. The component at 
the systemic velocity was detected in only the low- J tran- 
sitions (J = 5), implying that the absorbing gas producing 
it has a relatively low temperature. The blucshiftcd com- 
ponent, the strongest of the three overall, was detected in 
all of the transitions including those at high- J, which in- 
dicates a significant contribution from high-temperature 
gas. Its line profiles are asymmetric and could consist of 
two or more velocity components. The redshifted compo- 
nent, which is weaker than the others, was detected only in 
high- J transitions ( J^4), indicating that the temperature 
of the gas producing this component is higher than that 
producing the blucshiftcd component. IRAS 08572+3915 
is the first extragalactic object in which the fundamental 



band of CO has been resolved into individual lines, and 
where the lines themselves are velocity-resolved. 

5. Physical conditions of molecular clouds 

5.1. Temperature and column density 

5.1.1. Full coverage case 

In this and the following subsection we analyze the spec- 
trum in order to derive column densities and temperatures 
of the warm and cold components. In this subsection we 
make the simple assumption that the area covering fac- 
tor (C/) of the absorbing molecular gas is unity; in other 
words that for all transitions and at all velocities at which 
the CO absorbs, the absorbing gas completely covers the 
continuum source. 

In order to measure the intensities of the absorption 
lines, it is critical to know the continuum level accurately. 
Because the present spectrum is full of broad absorption 
lines, that level is difficult to be determined. We initially 
regard the local maxima of the continuum between each 
line as defining the continuum level. In that case the con- 
tinuum flux density decreases gradually with increasing 
wavelength at wavelengths longer than 5.00 {im. We fit 
this portion of the continuum with a linear function. The 
apparent continuum on the short-wavelength side of the 
spectrum shows significant structure and is inconsistent 
with a linear dependence on wavelength. We attribute this 
to two emission lines: HI Pf/3 and H2 v=0— >0 5(9) cen- 
tered at 4.65378 fim and 4.69462 ^m in the rest frame, re- 
spectively. Goldader et al (1995) and Higdon et al (2006) 
have reported the detection of H I and H2 emission lines 
at near and mid-infrared wavelengths from this galaxy. 
The lower and upper level energies associated with the 
H2 5(9) line are approximately 4,600 K and 7,200 K, re- 
spectively and hence it is not surprising that the line is 
observed in emission. We assume that each of H I Pf/3 and 
H2 5(9) lines have Gaussian profiles with the same veloc- 
ity widths as HI Br7 (241 km s" 1 ) and H 2 v=l->0 5(1) 
(430 km s" 1 ) observed by Goldader et al. (1995). The 
resulting best fit continuum spectrum synthesized from 
a straight line plus two Gaussian profiles is also plotted 
in Figure 2. Our analysis of the few CO lines coincident 
in wavelength with these emission lines assumes that the 
emission lines are behind the absorbing CO. This appears 
to us to be the most likely situation for the Pf/3 line, which 
arises in ionized gas. The H2 5(9) line occurs in molecular 
gas much hotter than the warm red- and blue-shifted gas 
producing CO absorption. It is possible, however, that it 
is physically associated with some of that gas and is in 
the foreground. However, the line is weak relative to the 
continuum and the continuum placement has only a small 
effect on the column densities of the one or two coinci- 
dent CO lines and no effect on the principal results of this 
paper. 

The equivalent widths of the CO absorption lines, mea- 
sured against the continuum level discussed above, were 
used to calculate the CO column density in the lower state 
of each transition. Since we found it impossible to defini- 
tively assign specific velocity intervals to the three velocity 
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components discussed earlier, only total line equivalent 
widths Wj across the entire profile were determined by 
numerical integration of the optical depth r, using 



Wj 



Io-h 
To 



X 2 



(1 



dv 



(1) 



where J is the rotational quantum number of the lower 
state, Iq is the continuum level, I\ is the intensity of the 
spectrum, and A and v is are wavelength and frequency of 
the line profile (but in the right portion of the equation 
A is the wavelength of the line center) . When stimulated 
emission is negligible compared to absorption, the column 
density in the lower state of a transition Nj is related to 
the equivalent width by 



N, 



1 



3 2 /A 2 W J > 



(2) 



in the optically thin case, where ire 2 jm e c 2 = 8.8523 x 
10 -13 cm, and / is the oscillator strength of the transi- 
tion. The oscillator strengths were taken from Goorvitch 
& Chackerian (1994). The measured equivalent widths 
and column densities are summarized in Table 1. 

On the assumption that the CO is in local thermody- 
namic equilibrium (LTE), the column density Nj in the 
rotational level J follows the Boltzmann distribution, 



N 



Z(T CX ) 



(2J + 1) exp - 



(3) 



where TV is the total column density, T cx is the excitation 
temperature, Z(T ex ) is the partition function which can 
be calculated by a polynomial interpolation provided in 
the HITRAN'08 database (Rothman et al. 2009), and 
Ej rj BqJ(J + 1) is the energy of the J level where B 
is the rotational constant: Bo/k = 2.765 K. It is useful 
to plot the population diagram, \a[Nj/(2 J + 1)] versus 
Ej/k. If LTE is valid, the data points should be fit by 
one or more straight lines whose slopes are 1/T ex - 

Figure 5 shows the population diagram for the observed 
CO lines (where N is determined over the full velocity ex- 
tent of each line). The diagram indicates that there are at 
least two components with different excitation tempera- 
tures. The excitation temperatures, determined by least- 
squares fits, are T e x,coid = 23 ± 1 K and T eXjWarm = 325 ± 5 
K. The CO column densities of these components calcu- 
lated from Equation 3 are 7V C o,coid = (5.74 ± 0.06) x 10 17 
cm" 2 and TVco.warm = (2.65±0.04) x 10 18 cm~ 2 for 23 
K and 325 K components, respectively. Assuming the 
standard cosmic abundance ratio Nco/Nn 2 ~ 1.8 x 10~ 4 
in dense molecular clouds (Dickman 1978), these values 
correspond to /V H2 , C oid = (3.19 ±0.03) x 10 21 cm" 2 and 
/VH 2 ,warm = (1.47 ± 0.02) x 10 22 cm" 2 , respectively The 
derived temperature of the warm component is far higher 
than that found toward Galactic interstellar medium lines 
of sight (Lutz et al. 1996), and similar to that observed to- 
ward the obscured ULIRG IRAS F00183-7111 by Spoon 
et al. (2004). 

5.1.2. Partial coverage case 

As discussed later in Section 5.4, the 5 (im continuum 
probably originates in dust at characteristic nuclear dis- 
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Fig. 5. Population diagram of CO absorption lines toward 
IRAS 08572+3915 NW (see text) for unity covering factor. 
Temperatures shown correspond to the slopes of the two 
straight lines. Error bars are statistical uncertainties (lcr), 
and do not include systematic errors. 

tances of a pc-scale. Krolik & Begelman (1988) concluded 
that when a dusty molecular torus of such dimensions is 
associated with an AGN the torus is composed of a large 
number of small high-density molecular clumps. The ge- 
ometry of the obscuring material in front of the nuclear 
continuum source in IRAS 08572+3915 is unclear, but 
may also consist of molecular clumps. Hence partial cov- 
erage of the continuum source should also be considered. 

Although in the previous subsection we derived the 
physical parameters of the gas clouds assuming that the 
covering factor is unity, that assumption may not be ap- 
propriate for the following reason. Generally speaking, 
analysis of the R- and P-branch transitions which origi- 
nate from the same J level should give the same column 
density of Nj. However, on the assumption of a cover- 
ing factor of unity, the column densities calculated for the 
R- branch at J=l and J=2 are systematically lower than 
those for the P-branch. On the other hand, the oscillator 
strengths of the P-branch are systematically larger than 
those of the P-branch. The most probable explanation for 
this mismatch is that the absorption lines of the P-branch 
(P(l), P(2)) are saturated but the absorbing gas does not 
fully cover the continuum emitting region, so that some 
unobscured continuum is observed; i.e., the covering fac- 
tor Cf in the line of sight is less than unity. 

The optical depths of the low- J lines can be estimated 
by the curve-of-growth method. In Figure 6, we show 
curves of growth, assuming for simplicity that an intrin- 
sic line profile can be represented by a Gaussian function. 
Three curves of growth are shown, for different velocity 
widths (FWHM = 2cr\/2hi2 = 300, 200, 120 km s" 1 ; 
where a is the dispersion of a Gaussian function) . Over- 
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plotted on the figure are the absorption-line data from 
Table 1. In the optically thin part of the curve the equiv- 
alent width is proportional to the column density, while 
the equivalent width saturates in the optically thick part. 
For each absorption line plotted, the vertical position is 
determined from the observed equivalent width. For each 
absorption line pair from the same J level, i.e., and 
P(l), R(2) and P(2), the horizontal positions have been 
shifted along the ir-axis, while keeping their ratio of /A 
values fixed, until the pair match a growth curve. Their 
final positions give the column density for rotational level 
J. The measured equivalent widths best match the model 
curve with FWHM = 200 km s _1 . The corresponding 
optical depth of the gaussian profile is calculated analyti- 
cally. 

Using this method the optical depths obtained for the 
J=l and ,7=2 levels of CO are ~ 7.8 and ~ 3.9, 
and Tfl( 2 ) ~ 2.6 and Tp( 2 ) ~ 1-7, respectively. These results 
suggest that the low- J lines are heavily saturated and thus 
that the covering factor (Cf) of the absorbing molecular 
clouds is closer to 0.6 than 1.0. Ideally, one should use the 
curve-of-growth analysis method for all absorption lines. 
However, this was not possible for J ^ 3 lines since their 
i?-branch lines were outside our wavelength coverage. We 
therefore applied Cf = 0.6 to all the lines, and derived the 
column density of each line by the corrected optical depth 

The equivalent widths and column densities of each transi- 



Optical Depth [t] 
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Fig. 6. Curves of growth assuming that an intrinsic line pro- 
file is a Gaussian function. Three lines are corresponding to 
the line full width half maxima (FWHMs) of 300, 200, and 
120 km s" 1 . A good match is obtained for FWHM = 200 
km s , and the peak optical depth in the case of FWHM = 
200 km s" 1 can be read off of the upper x-axis. 

tion line on the assumption of Cf = 0.6 are summarized in 
Tabic 1. For low- J lines, the column densities for Cf = 0.6 
are 2 — 4 times larger than those for Cf = 1.0, while for 
high- J lines, the column densities for Cf = 0.6 are slightly 
(< 2 times) larger than those for Cf = 1.0. 

Figure 7 shows the population diagram obtained on the 
assumption that Cf = 0.6; it should be compared with 
Figure 5. The derived temperatures and column densities 
for both cases are summarized in Table 2. For the cold 
component, the excitation temperature is T e x,coid = 24± 1 
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Fig. 7. Same as Figure 5, but assuming a covering factor of 
0.6. 

K, which is similar to the Cf = 1.0 case. The derived 
column density is iVco.coid = (1.98 ±0.04) x 10 18 cm" 2 , 
which is three times larger than that of the Cf = 1.0 
case. On the other hand, for the warm component the 
excitation temperature is T cx warm = 273 ± 2 K, which 
is slightly lower than that of the Cf — 1.0 case, and 
the column density is iVco.warm = (4.48 ± 0.04) x 10 18 
cm" 2 , which is slightly larger than that of the Cf = 1.0 
case. As a result, the corresponding H2 column densities 
are A^.coid = (1.10 + 0.01) x 10 22 cm" 2 and iV H2 ,™rm = 
(2.49 ± 0.02) x 10 22 cm" 2 for 24 K and 273 K compo- 
nents, respectively, when the standard value of [CO]/[H2] 
in dense clouds is adopted. 

We caution that because the CO absorbs over an ex- 
tremely wide velocity range the assumption that Cf is 
independent of velocity, based only on analysis of the low- 
J lines, whose equivalent widths have comparatively large 
contributions from the cold gas, is questionable. It also is 
not obvious that Cf =0.6 can be applied to all the lines. 
Thus there is some uncertainty in the estimates of the col- 
umn density. Note, however, that the difference between 
the results with Cf = 0.6 and 1.0 is not large for the warm 
component, which is the main topic of the present pa- 
per. Hence, we take Cf = 0.6 as the nominal value for the 
remainder of the paper. 

5.2. Origin of each velocity component 

Our observations have clearly shown the existence in 
front of the northwestern nucleus of IRAS 08572+3915 of 
three basic gaseous components associated with different 
velocity ranges centered at the systemic velocity and cen- 
tered at —160 km s" 1 and +100 km s" 1 relative to it, and 
with three widely different excitation temperatures, only 
two of which are apparent in Figure 7. In this subsection 
we discuss the origin of each velocity component. We have 



not been successful at cleanly isolating these components 
in velocity space for detailed analysis and comparison, us- 
ing standard deconvolution techniques, because the ve- 
locity profiles are too complicated. The velocity extents, 
velocity ranges, and degrees of overlap of the components 
vary with rotational level. Moreover, the profile of each 
velocity component is asymmetric, and each could consist 
of several sub-components. Nevertheless, the large tem- 
perature difference of the two dominant absorption com- 
ponents allows us to use Figure 4 and Figure 7 to crudely 
separate and simply characterize them. The much weaker 
third component differs sufficiently in temperature and 
velocity from the others that it can also be crudely char- 
acterized. 

The cold component at ^24 K revealed by Figure 7 
must correspond to CO absorption at velocities near the 
systemic velocity. As seen in Figure 4, this component was 
detected only in transitions from levels with J ^ 5, imply- 
ing a gas temperature around 30 K. We interpret this com- 
ponent as arising in interstellar clouds in the host galaxy 
located far from the nucleus, since its excitation tempera- 
ture is low and is similar to those of giant molecular clouds 
in the Galaxy. The presence of a strong 3.4 [ira carbona- 
ceous dust absorption feature toward IRAS 08572+3915 
(e.g., Imanishi, Dudley & Maloney 2006) suggests that 
a considerable fraction of the interstellar molecular gas 
along the line of sight is in diffuse clouds, as is the case 
toward the Galactic center (Whittet et al. 1997). In dif- 
fuse clouds almost all of the carbon (~ 99%) is atomic and 
there very little of it is in CO (~ 1%). Thus, if the gas 
traced by the observed cold CO component is entirely in 
diffuse clouds, the H2 column density could be two orders 
of magnitude larger than the previously estimated value. 
However, the present CO data do not in themselves al- 
low us to determine the origin of cold CO at the systemic 
velocity. Accordingly, the H2 column density of 1 x 10 22 
cm" 2 derived above is a lower limit. 

Similarly, the warm component at ~ 273 K, which is the 
dominant contributor to Figure 7, must be directly associ- 
ated with the strongest absorption component, centered at 
— 160 km s" 1 in Figure 4. This component was detected in 
all of the transitions and is strongest at J=3— 5, implying 
a gas temperature around 200 K. Outside of the Galactic 
center, the temperatures of molecular clouds do not exceed 
100 K (Dickman & Clemens 1983), except in very limited 
volumes such as those close to the luminous stars and 
shocked regions. It thus seems unlikely this large column 
density of warm negative velocity gas is located outside of 
the nucleus of IRAS 08572+3915 NW. We propose that 
this component is close to the nucleus and is heated by 
the central engine of the AGN. Meijerink & Spaans (2005) 
showed that molecular clouds exposed to X-ray radiation 
(X-ray dissociation regions) can be heated to tempera- 
tures such as those observed here (see also Maloney et 
al. 1996; Meijerink, Spaans & Israel 2006). The warm 
molecular gas may also be collisionally heated due to in- 
teractions of portions of the gas having different velocities. 
However, although such turbulent heating is a possibility, 
the negative velocities of this warm component indicate 
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Table 2. Summary of temperatures and column densities. 



Covering factor Component Temperature CO column density H 2 column density 

(K) (10 18 cm" 2 ) (10 22 cm' 2 ) 

C f = 1.0 Cold 23 ±1 0.574 ±0.006 0.319 ±0.003 

C f = 1.0 Warm 325 ±5 2.65 ±0.04 1.47 ±0.02 

C f =0.6 Cold 24 ±1 1.98 ±0.01 1.10 ±0.01 

C f = 0.6 Warm 273 ±2 4.48 ±0.04 2.49 ±0.02 



Notes. 

Col. (4): Total CO column density derived on the basis of the temperature in Col. (3). 

Col. (5): Total H2 column density, assuming the standard cosmic abundance ratio of Nco/Nn ~ 1.8 X 10" 4 (Dickman 1978). 



that the bulk gas motion is outward from the nucleus. 

The weak ±100 km s" 1 component seen in Figure 4 
docs not contribute significantly to Figure 7 because of 
its small column density. The equivalent widths of this 
component in the P(6) and P(16) line profiles are almost 
the same. From them we estimate that the temperature 
of the redshifted component is around 700 K and that the 
column density of CO is ~ 10 17 cm" 2 . If the source of the 
heating of this cloud is the central engine, then the red- 
shifted component would be expected to arise from clouds 
which lie closer to the central engine than those producing 
the blucshiftcd component. Outflows have been observed 
toward various types of AGN, but there are few examples 
of infall (e.g., Crenshaw et al. 2003; Muller Sanchez et 
al. 2009). We suggest that this velocity component traces 
the material fueling the central engine. 

5.3. Densities and thicknesses of the warm clouds 

The geometrical thickness of the absorbing cloud along 
the line of sight can be estimated from the column density 
derived above, if the gas particle density is known. As seen 
in the population diagram of Figure 7, the warm (~ 270 
K) component is well thcrmalizcd at least to the J=17 
level. This indicates that the population distribution of 
the CO molecules in the rotational states is controlled by 
collisions with other species, presumably mostly H2. 

A high gas density is required to thcrmalize CO up 
to the J=17 level. The density at which collisional de- 
excitation competes with radiative de-excitation (the crit- 
ical density; 7i. c (H2)) for CO J — > J — 1 is approximately 
n c (£l2) = 4 x 10 3 J 3 cm" 3 , as described in Kramer et al. 
(2004). The resulting value is n c(17 ^ 16 ) (H 2 ) = 2 x 10 7 
cm" 3 . Collisional transitions with A J > 1 arc much 
less likely than A J = 1 transitions and hence the criti- 
cal density is approximately the above value. This lower 
limit to the actual density is extremely high compared 
to the densities of typical molecular clouds in the Galaxy 
(n ~ 10 3 cm" 3 ). Such densities are, however, found in 
post-shock regions of molecular clouds, where molecular 
outflows from young stars collide with ambient molecular 
gas (e.g., Burton et al. 1988; Vannier et al. 2001), and 
may also occur where molecular clouds collide. 

Assuming that the absorbing cloud is homogenous, the 
thickness of the absorbing layer, Ad, is 



Ad< — ^- « 4x 10" 4 pc. (5) 
n c (H 2 ) 

The absorbing gas column thus appears to be extremely 
short along the line of sight. In general, dense molecu- 
lar clouds are expected be highly clumpy, as a natural 
result of the interplay between the energy generated by 
self-gravity of the gas and radiative cooling, dumpiness 
increases the size of molecular clouds with respect to the 
homogeneous case estimated above. It is rather difficult to 
estimate the volume filling factor, the fraction of absorbed 
volume in the beam area. Assuming that it is about 1- 
10%, similar to the values in highly clumped clouds in 
the Galaxy (Genzel et al. 1985), the thickness of the ab- 
sorbing molecular clouds would be 10-100 times greater; 
i.e., 0.004-0.04 pc. This is still much less than the thick- 
nesses of molecular tori assumed in previous studies of 
AGN (e.g., Jaffc et al. 2004; Wada & Tomisaka 2005). 

5.4- Extinction estimates and locations of emitters and 
absorbers 

As described in Section 5.1, we have detected the Pf/3 
hydrogen recombination line. The intensity ratios of hy- 
drogen recombination lines can be used to estimate the 
extinction by foreground dust, assuming "Mcnzel case -B" 
conditions (Osterbrock & Ferland 2005), for typical of H II 
regions. We thus compare the observed flux of Pf/3 (4.6538 
/Ltm, (3.2 ±0.5) x 10" 17 W/m 2 ), with that of Br 7 (2.1661 
/an, (0.7±0.2) x 10" 18 W/m 2 ) observed by Goldader et al. 
(1995). The flux ratio Pf/3/Br7 is ~ 46, whereas case B re- 
combination theory predicts a ratio of 0.58. The difference 
suggests a differential extinction of E(M — K) = 4.75 mag 
to the recombination line emitting region. According to 
the standard extinction curve model for Milky Way dust 
(Draine 2003), this value corresponds to Ay = 57 mag, im- 
plying a H column density of TVh = 1-07 x 10 23 cm -2 (using 
iV(H)/A v = 1-9 x 10 21 cm 2 mag" 1 ; Bohlin et al. 1978). 
The ratio could be considerably underestimated because 
the aperture used for our Pf/3 observation (0"60 x 0"46) 
is 100 times smaller than that for Br7 (3."0 x 9."0). Thus, 
the above column density is a lower limit. Case B condi- 
tions are found in H II regions, but not in dense ionized 
winds; we caution that the physical conditions in the gas 
producing the hydrogen recombination lines is unknown. 

From Section 5.1.2 the H 2 column density summed 
from all CO velocity components is Nn 2 ~4 x 10 22 cm" 2 , 
slightly lower than the lower limit estimated from the H 
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recombination lines. However, this estimate of iV(H2) is 
itself a lower limit in that, as discussed earlier, some or all 
of the cold CO may reside in diffuse gas. X-ray observa- 
tions toward IRAS 08572+3915 imply that the H column 
density to the AGN exceeds N K » 10 24 enr 2 (Risaliti et 
al. 2000). For comparison, the total hydrogen column 
density estimated from 9.7 /xm silicate dust absorption is 
N H ~ 1 X 10 23 cm -2 {t 9 . 7 /A v = 0.05-0.1; Roche k Aitken 
1984, 1985; r 9 . 7 = 4.2; Spoon et al. 2006), comparable to 
the lower limit to the H2 column density derived from the 
CO absorption. Wc conclude that the column density of 
hydrogen in front of the AGN, although highly uncertain, 
exceeds 10 23 cm -2 . 

Due to the large uncertainties in the above estimates, we 
cannot confidently assign relative locations along the line 
of site to the X-ray emitting region, recombination line 
emitting, infrared continuum emitting, and warm molec- 
ular absorbing regions based solely on estimates of fore- 
ground column densities of hydrogen. On the basis of 
energetics, however, one expects that the X-ray emitting 
gas must be interior to the recombination line emitting 
region. Both of those must be interior to the infrared 
continuum emitting regions at 5 and 10 microns, which 
themselves must be internal to the CO and silicate ab- 
sorptions, respectively. These assignments are in fact con- 
sistent with the column densities derived above if they 
are actual values rather than (in some cases) lower limits. 
Also these suggest that the infrared continuum emission 
comes, not from the central engine itself, but from a lo- 
cation somewhat further from the engine. One expects 
that the continuum emission in the M-band originates 
close to the inner edge of the dusty molecular clouds clos- 
est to the AGN. The dust closest to the AGN probably 
has a high temperature, roughly equal to the sublimation 
temperature around ~ 1500 K. This hot dust produces 
thermal emission whose spectral energy distribution peaks 
at near-infrared wavelengths, and also contributes signif- 
icantly at 5 (im continuum emission (Neugebauer et al. 
1987; Barvainis 1992; Kobayashi et al. 1993). The CO 
absorption lines would then be observed largely against 
continuum emission produced by this hottest dust near- 
est to the AGN, while the silicate absorption would occur 
largely against dust located somewhat more distant from 
it. 

6. Discussion 

Our spectrum toward IRAS 08572+3915 NW clearly 
separates the CO absorption into three velocity compo- 
nents; a cold component at the systemic velocity, a strong 
and warm component that is highly blueshifted (—160 
km s -1 ), and a weaker, but much warmer component that 
is highly redshifted (+100 km s -1 ). From the population 
diagram of the observed CO absorption lines, we deduce 
that the blueshifted absorbing gas is dense (> 10 7 cm" 3 ) 
and warm (~ 270 K) and has a CO column density of 
Nqo ~ 4.5 x 10 18 cm -2 , which is equivalent to a H2 col- 
umn density of N^ 2 ~ 2.5 x 10 22 cm~ 2 . The redshifted ab- 
sorbing gas appears to be produced by even warmer and 



presumably denser gas. These high temperatures and high 
densities, as well as the small column lengths (< 0.04 pc), 
indicate that these absorption components arc in close 
proximity to the central engine. 

The physical properties that we derive are similar to the 
recent results from near- to mid-infrared spectroscopy of 
other ULIRGs. Lahuis et al. (2007) detected strong ab- 
sorption bands of gaseous C2H2, HCN, and CO2 toward 
deeply obscured (U)LIRG nuclei. Their data suggest the 
presence of warm (200-700 K) and dense (> 10 7 cm -3 ) 
gas. Armus et al. (2007) argued that this gas occupies 
only a small fraction of the nuclear region (~ 0.01 pc) 
near the mid-infrared continuum source. However, be- 
cause all of their observations were made at much lower 
spectral resolution, they did not resolve individual veloc- 
ity components. Our observations at high spectral reso- 
lution shed additional light on the geometrical structure 
and kinematics of the molecular clouds in one ULIRG and 
our conclusions may apply to others as well. 

The extremely deep silicate absorption in some ULIRGs 
requires their nuclear sources to be deeply embedded in 
optically and geometrically thick material (Levenson et 
al. 2007). Such absorptions together with pointlike in- 
frared nuclear sources often are taken as signposts of the 
existence of AGNs (Imanishi, Dudley k Maloney 2006). 
However, Spoon et al. (2006) discovered crystalline sil- 
icate features toward IRAS 08572+3915 in addition to 
the deep amorphous silicate feature, which can be a 
sign of high-tempcraturc grain processing by the merger- 
triggered star-formation driven dust injection. 

Various kinds of theoretical models for the material lo- 
cal to the AGN and obscuring it have been proposed. In 
all of them the central engines are obscured by dusty, opti- 
cally and geometrically thick molecular clouds in toroidal 
structures; for example, a torus (Schartmann et al. 2005), 
a flared disk (Fritz et al. 2006), and discrete clouds 
(Elitzur k Shlosman 2006). Krolik k Begelman (1988) ar- 
gued that the obscuring torus is composed of high-density 
molecular clouds. Nenkova et al. (2002) concluded that 
the obscuring torus is composed of a large number of small 
but dense molecular clouds. Indeed, the thermal dust 
emission from the torus has been spatially resolved by 
recent IR interferometric techniques and shown to have a 
clumpy or filamentary dust structure (Jaffe et al. 2004; 
Tristram et al. 2007). Additional evidence for the exis- 
tence of the dust torus with many clumpy clouds comes 
from the detection of the silicate feature at ~10 (im in the 
spectral energy distribution of AGNs. The recent discov- 
eries of silicate emission from QSOs (Sicbcnmorgcn et al. 
2005; Hao et al. 2005), Seyferts (Hao et al. 2007), and 
LINERs (Sturm et al. 2005) by Spitzer provide strong ev- 
idence for a clumpy torus model. From a theoretical point 
of view, according to the 3-dimensional radiative transfer 
model of clumpy dust tori of AGN (Vollmer, Beckert k 
Duschl 2004; Beckert k Duschl 2004; Honig et al. 2006; 
Honig k Beckert 2007; Nenkova et al. 2008), the proper- 
ties of the individual self-gravitating clouds are speculated 
to be high density (n ~ 10 7 ~ 8 cm -3 ), large column den- 
sity (A H ~ 10 24 cm -2 ), and small volume filling factor (~ 
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0.03). Wada (2007) investigated the possible origins of 
molecular absorption, based on a 3-dimensional hydrody- 
namic model of the inhomogencous interstellar medium 
around the AGN. They demonstrated that the warm and 
dense molecular clouds, such as those discovered here via 
their CO absorption, can survive in the form of high- 
density clumps in the inhomogeneous interstellar matter 
around the AGN. 

Our high resolution spectroscopy has successfully de- 
tected warm dense gas with complex velocity structure, 
containing both highly blucshifted and highly rcdshiftcd 
components. Moreover, the CO spectrum suggests that 
the absorbing clouds have a large area covering factor 
(> 0.6) but a small volume filling factor (thickness < 0.04 
pc) . Because the velocity dispersions of the redshifted and 
blueshifted components are so large, it seems inconceiv- 
able that each exists in a single narrow sheet; hence each 
component must itself be composed of numerous well- 
separated thin sheets of dense gas that are detached from 
the continuum source. 

The redshifted and blueshifted components may trace 
the fueling of the central source and the ejection of ma- 
terial from it. However, most previous theoretical stud- 
ies (e.g., Wada & Tomisaka 2005) have concluded that 
the dominant velocity structure in the torus is system- 
atic rotation and turbulence. For IRAS 08572+3915 NW 
such models arc apparently inconsistent with our ob- 
served multiple, discrete velocity components. The pic- 
ture that the absorbing clouds in the torus rotate around 
the AGN cannot reproduce the observed CO absorptions 
in IRAS 08572+3915 NW for the following reasons: (1) 
Absorption by gas clouds in orbital motion around the 
AGN should reproduce a broad absorption profile cen- 
tered at the systemic velocity, because it will be make up 
of the integrated absorption of the entire line-of-sight ve- 
locity. The observed deep absorption lines with clear line 
separation into the blue-shifted and red-shifted compo- 
nents cannot be expected; (2) The redshifted absorption 
is extremely weak relative to the blueshifted absorption; 
and (3) The rcdshiftcd gas has much higher temperature 
than the blueshifted gas. 

We conclude that the velocity-shifted line components 
do not arise in a stable rotating torus-like gaseous struc- 
ture, but are more likely associated primarily with ejection 
and injection of a dense and warm gas close to the AGN. 
Our results suggest that, at least in this case, the environ- 
ment immediately surrounding the AGN does not appear 
to be as simple as that proposed in the unified scheme of 
AGNs. 

7. Summary 

We have observed and analyzed a high-resolution M- 
band spectrum of the obscured AGN, IRAS 08572+3915 
NW, which contains strong absorption lines of the fun- 
damental (y=l4— 0) band of CO. The main results and 
conclusions are as follows. 

(1) Absorption lines of 12 C 16 v=l<— are seen up to an 
excitation level of J=17. 



(2) The absorbing gas has three velocity components: 
a cold component at the systemic velocity; a strong 
and warm component that is highly blueshifted (—160 
km s _1 ); and a weaker, but much warmer component that 
is highly redshifted (+100 km s _1 ). 

(3) Assuming LTE in the absorbing molecular gas, the 
temperatures of the systemic and blueshifted components 
are 23 ± 1 K and 325 ± 5 K, and the CO gas column 
densities are Afco.cold = (5.74 ± 0.06) x 10 17 cm -2 and 
Aco,warm = (2.65 + 0.04) x 10 18 cm -2 , respectively, for a 
covering factor of unity. For dense molecular gas these 
values correspond to A^coid = (3.19 + 0.03) x 10 21 cm~ 2 
and A^.warm = (1.47 + 0.02) x 10 22 cm" 2 . In the case of 
a covering factor of 0.6, which seems more likely based 
on optical depth arguments, the temperatures of these 
components are 24 + 1 K and 273 + 2 K, and the CO col- 
umn densities are iVco.coid = (1.98 + 0.01) x 10 18 cm~ 2 
and AT C o,warm = (4.48 + 0.04) x 10 18 cm" 2 . These values 
correspond to 7V H2 .coid = (1.10 + 0.01) x 10 22 cm -2 and 
N H2 ,warm = (2.49 + 0.02) x 10 22 cm" 2 , respectively. The 
rcdshiftcd component has a higher temperature (~ 700 
K) than the blueshifted component and its column den- 
sity is considerably less (A*co,hot ~ 10 17 cm~ 2 , yielding 
A^H 2 .hot ~ 10 21 cm -2 ). The H2 column density associated 
with the cold gas at near systemic velocities may be con- 
siderably larger than the above values if the cold CO is 
largely in diffuse gas clouds. 

(4) From an estimate of the critical density, we demon- 
strate that the molecular gas with warm temperature 
must be very dense (n c (H.2) > 2 x 10 7 cm~ 3 ) and geo- 
metrically thin (< 0.04 pc) along the line of sight, even if 
it is highly clumpy. 

(5) The presence of discrete, highly blueshifted and highly 
rcdshiftcd absorption components close to the central en- 
gine docs not fit the standard unified model of AGNs. 
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